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ABSTRACT 

A massive black hole (MBH) in a galactic center drives a flow of stars into nearly radial orbits to 
replace those it destroyed. Stars whose orbits cross the event horizon rs or the tidal disruption radius 
rt are promptly destroyed in an orbital period P. Stars with orbital periapse r p slightly larger than 
the sink radius q = max(rs, r t ) may slowly spiral in due to dissipative interactions with the MBH, e.g. 
gravitational wave emission, tidal heating or accretion disk drag, with observable consequences and 
implications for the MBH growth rate. Unlike prompt destruction, the inspiral time is typically 3> P. 
This time is limited by the same scattering process that initially deflected the star into its eccentric 
orbit, since it can deflect it again to a wider orbit where dissipation is inefficient. The ratio between 
slow and prompt event rates is therefore much smaller than that implied by the ratio of cross-sections, 
~ Tp/q, and so only prompt disruption contributes significantly to the mass of the MBH. Conversely, 
most stars that scatter off the MBH survive the extreme tidal interaction ("tidal scattering"). We derive 
general expressions for the inspiral event rate and the mean number of inspiraling stars, and show that 
the survival probability of tidally scattered stars is *~ 1, and that the number of tidally heated stars 
( "squeezars" ) and gravity wave emitting stars in the Galactic Center is ~ 0.1-1. 

Subject headings: black hole physics — Galaxy: center — stellar dynamics — gravitational waves 



1. INTRODUCTION 

Evidence for the presence of massive black holes in the 
centers of most galaxies (e.g. Gebhardt et al. 2003), to- 
gether with recent and anticipated advances in observing 
capabilities, have focused interest on the observational im- 
plications of star-MBH interactions such as tidal disrup- 
' tion (Frank & Rees 1976; Frank 1978; Lightman & Shapiro 
1977; Syer & Ulmer 1999; Magorrian & Tremaine 1999), 
tidal scattering (Alexander & Livio 2001), gravitational 
' wave (GW) emission (e.g. Hils & Bender 1995; Sigurdsson 
, & Rees 1997; Freitag 2001, 2003), interaction with a mas- 
■ sive accretion disk surrounding the MBH (Ostriker 1983; 
Syer, Clarke & Rees 1991; Vilkoviskij & Czerny 2002), 
and tidal capture and tidal heating (Novikov, Petchik & 
Polnarev 1992; Alexander k Morris 2003, AM03). Much 
effort has been devoted to the study of tidal disruption, 
which plays an important role in the growth of low-mass 
MBHs (Murphy, Cohn & Durisen 1991; Freitag & Benz 
2002) and can provide a signature for the existence of MBH 
in galactic nuclei by the emission of tidal flares. 

Dynamical analyses indicate that most of the stars scat- 
tered into radial orbits originate at the MBH radius of in- 
fluence, r/j, where the enclosed stellar mass roughly equals 
the MBH mass m and the scattering is roughly isotropic. 
Event horizon crossing or tidal disruption is prompt; the 
stars plunge toward the MBH with a cross-section that 
scales as ~ r p (Hills 1975) and reach it in less than the 
initial orbital period, Pq, irrespective of orbital periapse, 
as long as r p < q. Thus, the star is destroyed in a short 
time, e.g. Po( r /i)~10 5 yr in the Galactic Center (GC). 

In contrast, the other processes listed above proceed 
gradually over an inspiral time to 3> Pq, as a small frac- 
tion of the orbital energy is dissipated every peri-passage. 
The inspiral time typically rises steeply with the periapse. 
In order for the extracted orbital energy to heat the disk, 
tidally heat the star or power a high luminosity of gravity 



waves, the star has first to decay into a short period orbit. 

Novikov et al. (1992) estimate that tidal capture by a 
MBH occurs for orbits with r p /r t <bo^3. It then follows 
that stars are scattered into tidal capture orbits at a rate 
bo— 1~2 times faster than that for prompt tidal disruption 
orbits. The orbital energy the star has to lose to circular- 
ize far exceeds its own binding energy, so it is likely that 
it will ultimately be disrupted (Rees 1988; AM03). This 
has led several authors (Frank & Rees 1976; Novikov et 
al. 1992; Magorrian & Tremaine 1999) to suggest that 
slow tidal inspiral may be at least as important as prompt 
disruption for feeding the MBH and for producing tidal 
flares. Simulations (Murphy et al. 1991; Freitag & Benz 
2002) indicate that prompt tidal disruptions supply be- 
tween ~0.15 to 0.65 of the total mass of a low mass MBH 
(m < 10 7 Mq) in a low-density galactic nuclear core. If 
the contribution of inspiraling stars were indeed as high as 
implied by the ratio of the cross-sections, this would have 
far-reaching implications: stars could supply most or even 
all of the MBH mass, thereby establishing a direct link 
between m and stellar dynamics on a scale of r^. 

However, a small initial periapse does not in itself guar- 
antee ultimate disruption. The star must also have enough 
time to complete its orbital decay. In this Letter we revisit 
the questions: what is the time available for orbital decay, 
and what is the inspiral event rate? 

2. CALCULATIONS 
2.1. Scattering into inspiraling orbits 

We follow the analysis of Miralda-Escude & Gould 
(2000, MG00) of the infall of a single mass population 
of stellar BHs into a MBH by diffusion into the loss-cone 
due to 2-body encounters in a Keplerian potential. 

Stars scattered from a volume dV to a periapse r p will 
spiral into the MBH in a time to(a,r p ), determined by 
the dissipation process. The scattering rate per volume of 
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stars into orbits with periapse <r p , dT(<r p )/dV, can be 
estimated for the steady state distribution function (DF) 
(MGOO Eq. 28). The total rate, T(< r p ) = J dV dT{< 
r p )/dV, is obtained by integrating over the volume be- 
tween an inner boundary ro where the stellar cusp is trun- 
cated (e.g by stellar collisions) to an outer boundary r c . 
The outer boundary is set by requiring that t <t p , where 
t p (a,r p ) = r p j (dr p /dt) is the orbit-averaged time for a 
change of order unity in the periapse by diffusion in veloc- 
ity space due to many small angle deflections. 

Typically, to rises with a while t p falls (Eq. 2), and 
so there is a critical semi-major axis a c (r p ) such that 
to{a Cl r p )/ t p (a c ,r p ) — 1. Stars originating from an orbit 
with a>o > a c do not have time, on average, to complete the 
inspiral. This simple condition is actually too restrictive 
because a(t) shrinks with time. Since stars are scattered 
off the orbit at a rate of ^t p 1 , the Poisson probability for 
avoiding this is w = exp(— s), where 
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The critical semi-major axis is obtained by solving s = 1 
for a c (r p ). Formally, T(< r p ) should be evaluated by 
the weighted integral J dV w dT(< r p )/dV over all space. 
However w falls off exponentially, and so r(< r p ) is well 
approximated by taking w — 1 and r c {r p ) = 3a c /2, the 
time-averaged radius on an e— >1 Kcplcrian orbit. 

We assume here a simple power- law stellar DF, n+ oc r~ a . 
In the limit e -> 1 (MGOO Eqs. 15-18, 21) 
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where Nh is the number of stars within r/,, Ai 
A(^p/^) 1/4 , A = m/M and 
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Note that r p /a= 1 — e = d 2 , where $ is the opening angle of 
the loss cone at r = a, so Eq. (2) is similar to the estimate 
^scattcr(a)^ 2 of Sigurdsson & Rees (1997). 

The volume contributing to T(<r p ) includes points at 
distance r far from the MBH, where to > t p for an orbit 
with periapse r p , but where to<t p is possible for r' p < r p 
due to higher dissipation at smaller periapse. We gener- 
alize T(< r p ) to account for this by solving s{r,r' p ) = 1 
(Eq. 1) for r' p (r), and setting r m (r) = mm[r p ,r' p (r)]. The 
scattering rate is then (Eqs. MGOO 15-18, 28), 

7 dr 
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where 7 = 7/2 — 2a, Ph = P{rh), r q = r c {q) is the maximal 
distance for prompt infall, and 

rF(3-a) 2 (10a- l)a\ 
2' 



B n = 



15 V 2 (a- 1/2)! 
Thus, the rate of successful inspiral events is 
T 4 (<r p ) =T m (<r p )-r m (<q). 



(a>-). (5) 
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The rate for prompt infall, T p (< r p ), is that at which 
stars are deflected into orbits with periapse < r p and reach 



r <r p at least once, but do not necessarily finish the inspi- 
ral. r p (<r p ) thus includes horizon crossing, tidal disrup- 
tion, inspiral and tidal scattering. The prompt infall time 
to ~ Pq does not depend on r p and so r m — r p at all r (for 
a<3; Eqs. 2, 15). The rate T p (<r p )^T m (<r p ) is then 

(7) 
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where Ei(x) = — J°° x dte 1 jt is the exponential integral. 

2.2. Inspiral and infall timescales 

Tidal heating. — We apply these results to "Hot 
Squeezars" (HS, AM03), tidally heated stars that dissipate 
the heat on the surface. We denote by a tilde dimensionlcss 
quantities in units of G=M i <=R+=\ (M*, i?* are the ini- 
tial stellar mass and radius). In these units r t — R{t)fh 1 / 3 '. 
The HS inspiral time in terms of b = r p /rt(0) is 

to - (2n) 2 ^m 2 ^b 6 P 1/3 /T(b 3 / 2 ) (m » 1), (8) 

where T is the tidal coupling coefficient. Two models are 
used to estimate the tidal heating: normal mode expan- 
sion (Press & Teukolsky 1977) and the affine stellar model 
(Carter & Luminet 1985). For the former, T is the leading 
multipole term, which is evaluated numerically for a so- 
lar model (Alexander & Kumar 2001). For the latter, we 
use the analytic approximation of Novikov et al. (1992). 
The tidal energy deposited in the star per peri-passage 
is AE = T/b 6 = const (AM03). The orbital evolution is 
(Eqs. 1-2, AM03 Eq. 4) 

(9) 
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s= to/ [2(3-a)t p (2o)] . 
The critical semi-major axis is 
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"Cold Squeezars" (CS, AM03) dissipate the tidal energy 
in their bulk and expand at constant effective temperature. 
The CS orbital evolution is calculated below numerically. 

Gravitational waves. — The time for inspiral by GW 
emission (b>q/r t ) is (Peters 1964) 

where c is the speed of light. Inspiraling GW emitters, like 
HSs, follow da/dtocy/a, so Eqs. (9)— (11) apply and 
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Prompt infall. — The mean infall time is t^ — Po/A, and 

a c .= r h . (15) 



logAiiV/, r h 
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3. RESULTS 

We apply these results to the GC by modeling it as a 
power-law DF with m = 2.6 x 10 6 M Q , r h = 1.8 pc and 
a = 1.8, based on the empirically derived mass model of 
Schodel et al. (2002). We assume a single mass popula- 
tion and Nh = rh. We consider 3 simple cases. (1) 1 M Q 
stars, a high mean mass motivated by theoretical argu- 
ments and observational evidence for a "top heavy" initial 
mass function in the GC (e.g. Morris 1993; Figer et al. 
1999). This model is used to test tidal inspiral and tidal 
scattering. (2) 0.1 M Q main-sequence (MS) stars, which 
are the most resilient against tidal disruption (Freitag & 
Bcnz 2002). This model is used to estimate GW inspi- 
ral from MS stars. (3) 0.6 M stars. Of these, 10% are 
white dwarfs (WD) with R+ — 0.01 Rq, as is expected in 
old, bulge-like stellar populations. The rest are MS stars, 
whose tidal disruption radius is too large for efficient GW 
emission. This model is used to estimate the rate of GW 
inspiral by WDs, for comparison with previous works. 

3.1. Survival probability of tidally scattered stars 

Tidal inspiral is complementary to tidal scattering 
(Alexander & Livio 2001), where stars narrowly escape 
tidal disruption by being scattered to wider orbits, af- 
ter suffering extreme tidal distortion, spin-up, mixing and 
mass-loss that may affect their evolution and appearance. 
Such stars eventually comprise a few percent of the pop- 
ulation within rh- We now show that their probability of 
survival is ~ 1 by comparing the inspiral and prompt infall 
rates. Table (1) lists the tidal inspiral rate Yi(<bor t ) (bo is 
the maximal possible pcriapse determined by r c (bor t ) =7"o) 
for 1 Mq HSs and CSs with the tidal heating estimated 
using cither normal mode expansion or the affine model 
(accounting numerically for the fact that t(R = b) <to for 
small b). The stellar cusp is truncated at r ^0.02 pc, the 
radius where stellar collisions destroy MS stars (Alexan- 
der 1999). The prompt disruption rate for this GC model 
is T p (< r t ) ~ 9x 10 _5 yr _1 , consistent with previous esti- 
mates, r p (<r t ) = 5xl0" 5 yr- 1 (Syer & Ulmcr 1999) and 
r p (<r t )~fewxlO~ 5 yr _1 (Alexander 1999). However, the 
squeezar inspiral (tidal capture) event rate is only ^0.05 
of the prompt disruption rate, and not < 2 times larger, as 
naively implied by the ratio of the cross-sections. 

The probability of successful inspiral for stars with pe- 
riapse between r t and r p is 

P '^ = rJT-i(<nr (16) 

while the tidal scattering survival probability is P s = 1— P». 
Our squeezar models have P s (< r t ) ~ 0.8 to P s (< b n r t ) ~ 
0.9. We thus confirm that P s ~l, as was anticipated from 
general arguments (Alexander & Livio 2001). The MBH's 
Brownian motion, neglected here, may further increase P s 
for loosely bound (6~&o) tidally scattered stars. 

A tidal scattering event is deemed "strong" if r p is 
within some adopted limit. Since P s ~l, the tidal scatter- 
ing and prompt disruption rates are related. The diffusive 
cross-section, modeled here, roughly scales as ^r 5 p , where 
£~(9-4a)/(8-2a) = 0.4for a = 1.8 (Eqs. 4, 15). Unbound 
stars with isotropic velocities have 8—1 (Hills 1975). Since 
most tidally scattered stars originate between the diffu- 
sive and the full loss cone (isotropic) regimes, where P^0 
(Lightman & Shapiro 1977), realistically 5~ 0.4-1. 



Table 1 

Inspiral in the Galactic Center 



Process™ 


ro 6 


r, 


n 


bo 


Li 


Pi 




pc 


y r ~ 1 








yr 


HS (S) 


2(-2) 


3(-6) 


0.2 


2.1 


170 L Q 


4(3) 


HS (A) 


2(-2) 


5(-6) 


0.4 


3.6 


150 L© 


6(3) 


CS (S) 


2(-2) 


4(-6) 


0.2 


2.8 


200 L© 


5(3) 


CS (A) 


2(-2) 


7(-6) 


0.4 


4.4 


170 L© 


7(3) 


MS GW 


6(-3) 


2(-7) 


0.2 


7r s 


2(35) erg/s 


1(1) 


WD GW 


4(-l) 


2(-7) 


0.04 


25 r s 


1(36) erg/s 


1(2) 



a S: Solar normal mode expansion. A: Affine stellar model. 
'The collisional destruction radius at tn = 10Gyr, defined 
here by e -'ff / 10t coii =0.1 (Murphy et al. 1991, Fig. 4a) 



3.2. Inspiral in the Galactic Center 

The observable implications of having on average n in- 
spiraling stars near the MBH can be estimated by consid- 
ering the properties of the leading (shortest period) star. 
The mean number is n(< feo^t) = J r ° rt dr p (to) (dTi/dr p ), 
where (to) <~ (2 + 2j)t (r p ,r c )/(3 + 2j) is the r-averaged 
value of to (tooc^/ao and Eq. 4). The mean inspiral time 
to = n(< bor t )/Ti(< bor t ) and the corresponding averaged 
Vao also define a and Po for "typical" HSs or GW emit- 
ters. For simplicity, we adopt this estimate also for CSs, 
although their to( a o) is different. The leading star com- 
pletes on average ti/to — n/ (n+ 1) of its inspiral time in 
N = n(n+2)to/2P orbits (AM03). Since a x =a /(n + l) 2 
and Pi = Po/(n + l) 3 , the total extracted orbital en- 
ergy is AEi = (rh/2ao)n (n + 2) and the tidal luminosity 
is Li (ti) - AEi/NPx = (m/a to)(n+ l) 3 . Typically, 
Li 3> i* (the intrinsic stellar luminosity) even for small n. 

Table (1) lists n, L\ and Pi for tidal and GW inspi- 
ral. We find that the GC contains on average ~ 0.1-1, 
squeezars, and that on average the leading squeezar is 
~ 5.5 m brighter than its normal bolomctric magnitude. 
The rate of WD inspiral derived here, Yi <~ 2 x 10 _7 yr _1 , 
agrees with the estimates of Sigurdsson & Rees (1997) and 
Freitag (2003), but the rate of GW inspiral by MS stars, 

<~ 2 x 10~ 7 yr _1 , is less than 0.1 of that estimated by 
Freitag (2003). The source of this discrepancy is unclear. 



4. DISCUSSION AND SUMMARY 

Different dissipation mechanisms may lead to the same 
outcome: orbital decay around a MBH in the presence of 
two-body perturbations. In this Letter we derive general 
expressions for estimating the inspiral event rate for any 
given inspiral time io( a o, r p) and orbital evolution a(t). In- 
spiral is a race between orbital energy extraction and two- 
body scattering. The probability that a star will reach the 
final, observationally interesting stage of a short period 
orbit is small unless it starts out on a tight enough or- 
bit. Since there are fewer stars close to the MBH, inspiral 
events are much rarer than prompt disruption events. 

We applied these results to the GC, using simple single- 
mass models to represent the stellar cluster. These repro- 
duce the prompt disruption and the WD inspiral event 
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rates that were independently estimated by previous stud- 
ies. We find that (1) The survival probability of tidally 
scattered stars is 0.8-0.9. (2) The rate of tidal scat- 
tering scales with the prompt tidal disruption rate as 
T p(< r t)[{r p /r t ) s - 1] with 5~0.4-l. (3) The contribution 
of slow tidal inspiral in the GC to the total tidal disrup- 
tion rate is only <~5%, and not 100-200% as proposed by 
previous studies. (4) The GC contains on average ~ 0.1-1 
squcczars at any given time, with a tidal bolometric lumi- 
nosity > 150i* on <~ 5 x 10 3 yr orbits. (5) There are on 
average ~0.1 GW emitters in the GC at any given time. 

The uncertainties in these estimates can be addressed 
by more detailed modeling of the stellar cluster, which 
should include a realistic, multi-mass stellar DF and take 



into account mass segregation. Other dynamical processes 
not considered here, such as resonant scattering (Rauch 
& Tremaine 1996), deviations from spherical symmetry 
(Magorrian & Tremaine 1999), chaotic orbits in triaxial 
systems (Poon & Merritt 2002), or the effects of massive 
pcrturbcrs (Zhao et al. 2002) may increase the inspiral 
event rates. The MBH's Brownian motion will not have a 
large effect, as typical inspiral orbits originate well within 
r h , where the stars follow the MBH (Reid et al. 2003). 

We thank M. Freitag and M. Morris for helpful discus- 
sions. TA is supported by ISF grant 295/02-1, Minerva 
grant 8484, and a grant by Sir H. Djangoly, CBE, of Lon- 
don, UK. 
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